Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
海老原 健一
JAERI-Research 2005-004, 121 Pages, 2005/03
JAERI-Research-2005-004.pdf:19.79MB
本報告書は、まず、格子ガス法及び格子ボルツマン法の二相流体モデルの二相流シミュレーションへの適用妥当性及び有用性を検討した。その結果に基づき、格子ボルツマン法の二相流体モデルであるHCZモデルを、基本的かつ重要な二相流流動様式の1つである水平層状二相流のシミュレーションへ適用した。その結果、HCZモデルの界面は、Kelvin-Helmholtz不安定性を満足することが確認された。またTaitelとDuklerの理論的流動様式線図(T-D線図)の界面成長に関する曲線が再現された。さらに、幅が狭い流路では、流れの三次元性の影響が顕著となり、T-D線図が与える以上の低密度相のみかけ流速が界面成長に必要であることが明らかとなった。次に、より複雑な界面現象を伴う液滴発生のシミュレーションを行い、Ishii-Grolmesの実験相関式を、実験データの分布範囲において、再現できることがわかった。
海老原 健一
筑波大学大学院システム情報工学研究科博士学位論文, 134 Pages, 2004/09
本論文では、格子法の気液モデルの有用性を検討し、さらに格子ボルツマン法気液モデル(HCZモデル)を、格子ボルツマン法の適用例が見られていない水平層状二相流の三次元シミュレーションへ適用した。その結果、以下のことが明らかになった。(1)HCZモデルによってシミュレーションされる二相界面が、Kelvin-Helmholtz不安定性理論を、二次元及び三次元の場合について満足することがわかった。(2)矩形断面流路中の界面成長シミュレーションにおいて、界面成長と流動状態の関係が、理論的考察によって得られるTaitelとDuklerの流動様式線図を再現できることがわかった。(3)流路幅が異なる場合の界面成長シミュレーションによって、流路幅が狭くなった時、流れの三次元性の影響が顕著となり、界面成長に必要な低密度相の流量が、理論的流動様式線図から得られる値より大きくなることが明らかとなった。(4)界面成長より複雑な界面現象を示す液滴発生シミュレーションにおいて、液滴発生とその時の流動状態の関係が、IshiiとGrolmesによって提案された実験相関式を、ほぼ再現できることがわかった。
海老原 健一; 渡辺 正
ながれ, 23(4), p.253 - 261, 2004/08
一成分二相格子ボルツマン法を水平層状二相流における液滴発生のシミュレーションに適用した。幾つかの流動状態において、液滴が界面上からちぎれる様子が観察された。また同時に、各々の場合において、流動状態を特徴付ける無次元数を測定した。シミュレーションに対する計算格子の長さや幅の影響を調べた後、液滴発生と測定された無次元数との関係を、IshiiとGrolmesによって実験結果に基づいて提案された相関式と比較した。その結果、シミュレーション結果は、実験データと相関式と差の範囲において、相関式を再現していることがわかった。
海老原 健一; 渡辺 正
日本機械学会論文集,B, 70(694), p.1393 - 1399, 2004/06
一成分二相格子ボルツマン法によって、その幅が高さより狭い矩形断面流路内における水平層状二相流のシミュレーションを行った。3種類の異なる流路幅の場合について、二相間の界面成長を測定した。さらに二相流動を特徴付ける無次元数を測定し、それらと界面成長との関係とTaitelとDuklerによって提案された流動様式線図とを比較した。流路幅が狭くなるに従い、シミュレーションにより得られた界面の成長・非成長の境界が大きくなることがわかった。
柴本 泰照; 久木田 豊*; 与能本 泰介; 安濃田 良成
Proceedings of 2nd Japan-Korea Symposium on Nuclear Thermal Hydraulics and Safety (NTHAS-2), p.374 - 380, 2000/00
原研で実施されたROSA/AP600実験において、静的余熱除去系作動後にコールドレグ内で最大温度差約150Kの顕著な温度成層が観察された。本報では、大気圧装置の実験を行い、このような成層形成時(流動冷水に対向して侵入する熱水の挙動)についての実験及び解析結果を報告する。実験では、実機相当の密度差を模擬するために、作動流体に飽和塩水と淡水を使用した。解析では、一次元運動量バランスを仮定したモデルを提案し、実験データ及び文献値との比較によりその妥当性を検証した。その結果、侵入速度が遅い領域では壁面及び密度界面の摩擦効果を考慮する必要があることを示した。またトタールの侵入流量については、自由界面の不安定性発生条件により説明できることを示した。
中村 秀夫; 近藤 昌也; 伊藤 高啓*; 辻 義之*; 久木田 豊*
Bulletin of the Research Laboratory for Nuclear Reactors, (SPECIAL ISSUE 2), p.35 - 46, 1997/00
超音波流速分布計測器(UVP)は、トランスデューサから超音波パルスを周期的に発振し、液相に混入したトレーサからの反射波とのドプラー信号から、液相の流速分布を計測する。本UVPを用いて、水平開水路の波状流での界面波形と液相流速分布の同時計測、並びに円管内に設置した小円筒後流に生じる渦列の特性の計測を試みた。水平波状流では、超音波が水面で全反射して界面・流路底間での多重反射を生じ、計測が困難であったが、超音波吸収体の設置等の工夫で、2次元的流速分布と水面形状の同時計測に成功した。円管内小円筒の後流に生じる渦列の計測では、低レイノルズ流で、渦周期の計測と共に、渦強度の減衰挙動を平均流速分布から捉えることに成功した。これらの計測結果をまとめる。
伊藤 和宏*; 辻 義之*; 玉置 昌義*; 中村 秀夫; 近藤 昌也; 久木田 豊
日本原子力学会誌, 39(8), p.669 - 680, 1997/00
被引用回数：0 パーセンタイル：0.01(Nuclear Science & Technology)水平水-空気波状流における気液界面せん断応力を実験的に評価した。実験では、高さ0.7m、長さ12mの矩形ダクトにおいて水の見かけ流速を0.3m/s一定とし、空気の見かけ流速を4.2~6.8m/sに変化させた。水面にh界面波が生じ、その振幅は空気流量又は入口からの距離と共に増加した。プレストン管法を用いて測定した気相部壁面摩擦係数は、この界面波の影響によって滑面に対するBlasius式の予測値と異なったが、その差は
30%であった。気液界面せん断応力は、壁面せん断応力、圧力勾配、水位及び水位勾配の各計測値を気相の運動量式に代入することで評価した。その結果、気液界面せん断応力は界面波の振幅と共に増加し、気液界面摩擦係数のBlasius式の予測値との比は4~10であった。この比は空気の見かけ流速と入口からの距離に比例して増加した。
浅香 英明; 久木田 豊
Journal of Nuclear Science and Technology, 33(9), p.696 - 702, 1996/09
被引用回数：2 パーセンタイル：24.74(Nuclear Science & Technology)加圧水型原子炉(PWR)の小破断冷却材喪失事故(LOCA)時におけるホットレグ内の流れは、重力により水と蒸気が分離した水平層状二相流で特徴づけられる。ここではホットレグ水位は、蒸気発生器(SG)内ないしは加圧器内への蒸気による水のキャリオーバー量に影響する重要なパラメータである。小破断LOCA時のホットレグ内層状二相流は、水流量が増加するにしたがい常流から射流に遷移し、これに起因する典型的な水位低下の生じることがLSTF実験により示された。Gardnerは、ホットレグ入口前後で各相のエネルギーが保存されると仮定し、常流から射流への遷移モデルを開発した。しかし、このモデルによる予測結果とLSTF実験結果との一致は、定量的に不十分であった。筆者らは、Gardnerのモデルにホットレグ入口の形状エネルギー損失項をとりいれた改良モデルを開発した。この改良モデルは、LSTF実験で見られたホットレグ内の常流から射流への遷移条件を良く予測した。さらに射流条件下の水位は、ホットレグ内における運動量保存則で予測できることを明らかにした。
中村 秀夫; 久木田 豊; 田坂 完二*
Journal of Nuclear Science and Technology, 32(9), p.868 - 879, 1995/09
被引用回数：4 パーセンタイル：43.23(Nuclear Science & Technology)LOCA時の1次系内冷却材分布は、水平配管内の流動様式に依存する。特に、波状流からスラグ流への遷移は、界面波が気送流れを塞ぐ液スラグに成長する結果生じ、界面摩擦や流路の圧力損失の増加を伴うため、その予測は重要である。ところが、スラグ流遷移の予測には、波状流での水位予測が必要であるものの、予測に必要な気液界面摩擦係数は、スラグ流遷移時のような液深の深い波状流(ボイド率
~0.6)でのデータや予測式が大変少なく、壁摩擦係数などで代用されることが多い。そこで、TPTF装置の水平円管テスト部(直径87及び180mm)を用いた、高圧(3-9MPa)の水/蒸気二相流実験を行い、スラグ流遷移境界近傍の波状流について調べ、界面摩擦係数を与える実験式を求めた。実験式はKelvin-Helmholtzの界面不安定性理論に基づくパラメータで構成され、TPTF実験結果での配管口径及び圧力依存性をよく表すことができる。
浅香 英明; 久木田 豊; 斉藤 誠司*
The 3rd JSME/ASME Joint Int. Conf. on Nuclear Engineering (ICONE),Vol. 1, 0, p.131 - 136, 1995/00
加圧水型原子炉(PWR)の小破断冷却材事故(LOCA)時におけるホットレグ内の流れは、重力により水と蒸気が分離した水平層状二相流で特徴づけられる。ここでホットレグ水位は、蒸気発生器(SG)内ないしは加圧器内への蒸気による水のキャリオーバー量に影響する重要なパラメータである。小破断LOCA時のホットレグ内層状二相流は、水流量が増加するにしたがい常流から射流に遷移し、これに起因する典型的な水位低下の生じることがLSTF実験により示された。Gardnerは、ホットレグ入口前後で各相のエネルギーが保存されると仮定し、常流から射流への遷移モデルを開発した。しかし、このモデルによる予測結果とLSTF実験結果との一致は、定量的に不十分であった。筆者らは、Gardnerのモデルにホットレグ入口の形状エネルギー損失項をとりいれた改良モデルを開発した。この改良モデルは、LSTF実験で見られたホットレグ内の常流から射流への遷移条件を良く予測した。さらに射流発生時の水位予測及び気液の相間摩擦が水位変化に及ぼす影響についても論じられている。
川橋 正昭*; 細井 健司*; 平原 裕行*; 福井 大俊*; 文沢 元雄
The 4th Triennial Int. Symp. on Fluid Control,Fluid Measurement,and Visualization; FLUCOME 94, 0, p.691 - 695, 1994/00
高温ガス炉空気浸入事故時の流動挙動解明の一環として、本研究では小口径開口部を通るヘリウムと空気の置換流について可視化観察を行った。可視化はレーザライトシート法を用いた。その結果、水平開口において捩れ構造の対向置換流が生じること、傾斜開口においては分離した二相流構造となることが明らかとなった。またこれらの可視化観察により、置換流量と傾斜角との関係を考察した。
大貫 晃; 秋本 肇; 村尾 良夫
Journal of Nuclear Science and Technology, 29(3), p.223 - 232, 1992/03
PWRホットレグでの気液対向流制限に対する二流体モデルコードのための界面せん断力モデルを、大貫らによる対向層状流での界面せん断力モデルをベースに開発した。まず、定常の包絡線モデルにより検証された大貫らのモデルが二流体モデルコードによる動的な計算において有効に機能するか否かを、代表的な二流体モデルコードであるTRAC-PF1/MOD1を使い解析した。その結果、大貫らのモデルはホットレグ内スラグ流での適切な界面せん断力モデルと組み合わせる必要のあることがわかった。ホットレグベンド部でのモデル及びホットレグ水平管内でのまき波領域に対するモデルをスラグ流モデルとして提案した。種々のスケール、圧力及び流体の種類(内径:0.025-0.75m、圧力:0.1-7.1MPa、空気-水または蒸気-水)のもとでの実験データにより本モデルを検証し、本モデルの有効性を確認した。
与能本 泰介; 田坂 完二*
Int. J. Multiph. Flow, 17(6), p.745 - 765, 1991/00
被引用回数：50 パーセンタイル：87.39(Mechanics)成層二相流領域から小破断口への液体及び気体のエントレインメントについて理論的及び実験的な検討を行なった。すべての破断形状について簡単に使用する事ができるように、以前に得られた理論式を破断流量、破断クオリティ及び水位の関係を示す式に修正した。さらに修正された式を常温大気圧条件で行なわれた実験データで評価した。実験結果は破断口周辺の流動条件が対称的である場合、本修正式とよく一致した。実験で観測された渦や横流れや波状流の影響については、本モデルを基にしてその相関関係を明確にした。この相関式を用いる事により高圧の水蒸気実験を含む様々な文献のデータを予測する事ができた。
片山 二郎; 中村 秀夫; 久木田 豊
Proc. of the Int. Conf. on Multiphase Flows 91-TSUKUBA,Vol. 1, p.7 - 10, 1991/00
LSTFを用いた小破断LOCA実験において、ホットレグ中に生じる層状流の液相流速が気液界面に生じる波の伝播速度を上回る、いわゆる射流が生じる場合があることが観察されている。本研究では、蒸気発生器入口配管を模擬した上昇管をもつ小口径水平配管を用いて大気圧下の水・空気層状流について射流の生じた場合の流動様式遷移に関する実験を行った。実験では、水平管内あるいは上昇管入口において跳水が生じ、跳水の下流側は常流となり層状流から波状流あるいはスラグ流へ遷移した。特に高い水流量では、射流から直接スラグ流へ遷移し、明瞭な跳水は認められなかった。この実験から以下の点が明らかになった。(1)水平管内でスラグ流への遷移が生ずる液流量は、常流と射流の場合で異なる。(2)しかし、スラグ流へ遷移する直前の条件での局所ボイド率と気液間相対速度は、射流・常流のいずれの場合も同一のモデルにより予測できる。
久木田 豊; 浅香 英明; 三村 裕一*; 安濃田 良成; 石黒 美佐子; 根本 俊行; 田坂 完二
JAERI-M 90-053, 22 Pages, 1990/03
現在開発の途上にあるRELAP5/MOD3コードの性能評価のため、ROSA-IV/TPTF装置による高圧(7MPa)、大口径(0.18m)水-蒸気水平二相流実験を行った。試験部内のボイド率に関する解析結果と実験結果の一致は、既存のRELAP5コード(MOD2)にくらべ著しく改善された。これは、層状流の発生限界に関するモデルならびに相間摩擦の計算方法がMOD2にくらべ改良された結果である。
大貫 晃
JAERI-M 85-219, 19 Pages, 1986/01
環状流および層状流の流動様式を示す円管内での空気-水対向流制限(CCFL)のデータを使い、TRAC-PF1/MOD1の予測性能の評価を行った。形状パラメータとして長さ、内径、傾斜角度をとり、加圧水型炉のホットレグを模擬した傾斜管付き水平管のデータに対しても評価を行った。環状流の相関摩擦係数に対しては、Wallis型およびBharathan型の2種類の相関式を調べた。データとの比較から、Wallis型の相関式のほうが落水が止まる点(バイパス点)をよく予測した。一方落水を生じる領域ではBharathan型のほうがよりよい予測を示すものの、正確な予測のためには円管の上下端に付加的な形状損失係数が必要なことがわかった。層状流およびホットレグ模擬管のデータに対しては、バイパス点を与える空気流速を過小評価するものの落水を生じる領域に対してはよい予測を与えた。
